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AB3TRA CT 

Empirical Orthogonal Function (EOF) representation of the 
synoptic forcing is combined with past meridional and zonal 
displacement s (persistence) to forecast tropical cyclone 
recurvature at 36 , 54 and 72 h. Recurvature is defined 

following Leftwich (1978, 1982) : Recurvature is a net 

displacement northward of 315° luring the forecast period or 
the attaining of northeastward motion during the 12 hours 
prior no the end of the forecast period. The combination of 
persistence and EOF coefficients consistently out-performed 
the individual methods for forecasting recurvature in both 
the dependent and independent data samples. Evaluation 
scores were the Brier p-score, the Heidke skill score and 
percent correctly forecast. 

Seven forecast aids used by Joint Typaoon Warning Center 
(JTWC) Guam, were tested for ability to forecast tropical 
cyclone recurvature at 24, 48 and 72 h. Using the skill 

scores described above for the 1979-31 typhoon seasons, 
performance indicates large annual variations of skill 
between the models. 



4 



•' v * t 





TA3LE OF CONTENTS 



I. IN TRODUCTION 10 

II. EOF ANALYSIS OF SYNOPTIC FORCING APPROACH TO 

RE CURVATURE 13 

III. DEFINITION OF RECURVATURE 23 

IV. REGRESSION ANALYSIS 30 

A. 36 H REGRESSION CASE ANALYSIS 38 

3. 54 H REGRESSION CASE ANALYSIS 41 

C. 72 H REGRESSION CASE ANALYSIS 44 

V. TESTING THE REGRESSION ANALYSIS . 47 

VI. CONCLUSIONS AND JUSTIFICATION FOR CONTINUED 

RESEARCH . 59 

APPENDIX A. MODEL COMPARISONS 62 

A. PURPOSE 62 

B. THE DATA 62 

C. MODEL DESCRIPTIONS COMPARISON TESTS 63 

D. SCREENING PROCEDURE FOR MODEL TESTING 66 

E. SUMMARY 71 

LIST OF REFERENCES 73 

INITIAL DISTRIBUTION LIST 76 



5 



LIST OF FI GORES 



Figure 1. Grid Arrangement Used in Statistical- 

Climatological Track prediction 13 

Figure 2. Eigenvectors 1-9 Multiplied by 100. The X 
Indicates the Location of the Tropical 
Cyclone, Grid Point 70 19 

Figure 3. Eigenvalues of the 500 mb D-value Fields 
(circles) and the Monte Carlo Generated 
Eigenvalues (triangles) Based on 100 
Simulations of Random Data 22 

Figure 4 . Examples of Idealized Recurvature (1,2), Non- 
recurvature (3) and a Not Considered Case (4) 24 

Figure 5. Typhoon Trix Case Illustrating an Apparent 

Recurvatura Case at 36 h, and Non-recurvature 
at 5 4 and 72 h 25 

Figure 6. Typhoon Polly Case Illustrating Recurvature 

at 36 and 72 h and Non-recurvature at 54 h . 25 

Figure 7. Typhoon Cera Case Illustrating a Common 

Recurvature Case 26 

Figure 8. Typhoon Alice Case Illustrating a Non- 
recurvature Case 26 

Figure 9. Typhoon Anita Case Illustrating the Exclusion 
or the 1 2 GMT 22 July to 00GMT 23 Julv Periods 
Due to the Eastward Motion 27 

Figure 10. Probability of Recurvature (504 Case Sample 

Climatology) 29 

Figure 11. Storm Tracks Relative to the Initial Position 
for All Storms with an Orthogonal Coefficient 
1 that is: (a) 3etween -9 and -30: or (b) 

Between +9 and +30. Each Cross Indicates a 
12-h Displacement (Shaffer, 1982) 37 

Figure 12. Heidke Skill Scores with the Value of the 
Predictand from the PSS Equations that is 
Defined t c be Recurvature Varied from 0.05 
tc 0.95. 54 



O 



LIST OF T ABLES 



TABLE 

TABLE 

TABLE 

TABLE 

TA3LE 

TABLE 

TABLE 

TABLE 

TABLE 

TABLE 

TABLE 

TA3LE 

TA3LE 



I. Eigenvalues and Cumulative Percent Explained 

Variance cf the Normalized D-Value Fields at 
500 mb (Shaffer, 198 2) 21 



II. Number of Recurving (rev). Non- recurving 

(nrev) , No History (nohist), Not Considered 
(notcon) and -he Probability of Recurvature 
(prob) for the 504 Tropical Storm Sample . . 2g 



III. Potential Predictors Used to Develop the 
Regression Equations for Each of the Three 
Time Intervals 32 



IV. Regression Equations for 36 h Recurvature 
for Each cf the Three Methods Tested 
(Dependent Sample), and Mean Values and the 
standard Deviation from the Means (STD) . . 



34 



Regression Equations for 54 h Recurvature as 
Described in Table IV 35 



VI. Regression Equations for 72 h Recurvature as 

Described in Table IV 35 



VII. Summary Cf SteDwise Linear Regression 

Results for 36 h (25 9 cases) R Squared (RSQ) , 
Increase in R Squared (RSQ Incr) and F to 
Enter Or Remove 41 



VIII. Summary of Stepwise 
Results for 54 h (20 6 
Table VII 



IX. Summary of Stepwise 
Results for 72 h (138 
Table VII 



Linear Regression 
cases) as Described in 
43 



Linear Regression 
cases) as Described in 
44 



XI. Predicta nd Statistics for the 36, 54 and 72 
h Before Filtering. The Number of Cases, 

Mean, Maximum and Minimum Values (max, mm) 
and the Variance of the Predicted Values for 
the Dependent and Independent Data Sets ... 51 



XII. Predictar.d Statistics for 36, 54 and 72 h 

After Filtering as Described m Table XI. . . 52 



XIII. Heidke Skill Score Using the Sample 

Ciimatolcqy as a Standard of ComDarison for 
Defining Recurvature Versus Noa-Recurvature 
for the Filtered Predictands (as in Table II) 

(A HSS of 1.0 is Perfect) 52 



XIV. Predicta nd Distrioutior. and Cumulative 

Percentage Distribution for 36 h. Cumulative 
Percentage is in Parens 



55 



TABLE XV. 



Predictand Distribution and Cumulative 
Percentage Distribution for 54 h. Cumulative 

Percentage is in Parens 56 

TABLE XVI. Predictand Distribution and Cumulative 

Percentage Distribution for 72 h. Cumulative 
Percentage is in Parens 56 



TABLE XVII. P-Score Statistics for 36, 54 and 72 h 

Dependent and Independent Data Sets After 
Filtering of the Predictand. (A Score of 
0.0 is Perfect) 57 



TA3LE XVIII. Percent Correctly Forecast Events frcm the 
Regression Method and Defining Recurvature as 
a Predictand Exceeding the Climatology 
Probabilities of Recurvature. (A Score of 
100 is Perfect) 58 



TABLE XIX. Number of Recurving (rev), Non-recurving 

(nrcv) , and Probability of Recurvature (prob) 
ror Tropical Cyclone Seasons 1979, 1980 and 

1981. . . . I 67 



TABLE XX. Heidke Skill Score for Models for 1979, 

1980, and 1981 {A Score Of 1.0 is Perfect, 

and a (*) Indicates too Few Cases to 

be Ranked) 69 

TABLE XXI. Erier ? - Score Statistics For 24, 48 and 72 
h Models for 1979, 1 980 and 1981 (A Score of 
0.0 is Perfect, and a (*) Indicates too Few 
Cases to be Ranked) 70 

TABLE XXII. Percent Correct Forecast Statistics for 

Models for 1979, 1980 and 1981 (A Score of 
100 is Perfect, and a (*) Indicates too Few 
Cases to be Ranked) 71 



3 



AC KN0WLED3 EM ENT 



I wish to thank: 

• Professor Russ Elsberry for his guidance, encouragement 
and patience given throughout the course of this 
research. He was there to assist <than assistance was 
needed and has taught me much moce than weeds can 
express. Thank you very much; 

• Professor R. L. Haney for his review of the thesis; 

• Dr. Ted Tsui of the Naval Environmental Prediction 
Research Facility for Droviding the data analyzed in 
the appendix of this thesis ; 



The people of the W 
their patience; 



R. Church computer center for 



Al Shaffer for his explanations, outdance, friendshiD 
and most imoortant, takincr the time to sit and discuss 
the mar.v problems and questions that were raised during 
the past nine months and; 



A very special thanks to Mary Rose, Abby and Tommy, for 
their understanding, encouragement and prayers during the 



past two years. 



9 



I. 



INTRODUCTION 



Tropical cyclone recurvature is one of the most 
difficult events to forecast for even the most experienced 
forecaster. The ability to forecast recurvature correctly 
could annually save both the military and civillian 
communities millions cf dollars, thousands of man hours and 
the lives of many victims. There are few reliable rules 
that the tropical cyclone forecaster can follow for 
forecasting recurvature cr non -rec urvature in either the 
near and the extended time frames. Riehl and Shafer (1944) 
were among the first to summarize guidelines for forecasting 
recuivature cf Atlantic tropical cyclones. Riehl and Shafer 
defined recurvature as "any major alteration of a storm 
track" which is quite non-descript. Riehl and Shafer listed 
several factors affecting tropical cyclone motion, including 
the effects of the upper-level winds and the interaction 
with troughs in the westerly flow to the west of the 
tropical cyclone. They also described various synoptic 
situations to assist the forecaster in predicting the storm 
path. The Joint Typhoon Warning Center { JTWC) Guam defines 
recurvature as 
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"The turning of a tropical cyclone from an initial path 
toward the west and northwest xo the north then northeasx" 

which is also quite non-descript. 

More recent studies by George and Gray (1977), Char., et 
al (1980) and Shaffer (1982) have shown that the tropical 
cyclone motion is largely dependent upon the surrounding 
steering flow represented by hei ght/temperature fields. 

The purpose of this study is to use a simple 

representation cf the synoptic forcing by Empirical 
Orthogonal Functions (EOF's) plus persistence-related 
variables to predict tropical cyclone recurvature or 
non-recurvature. The time increments of 36, 54 and 72 h 
for forecasting recurvature were selected for several 
reasons. First, a minimum time period of 36 h was thought 
to be of greatest value and to provide the most confident 
forecast for the shore establishment, air wing squadrons and 
ships presently in port. Secondly, The upper-air soundings 
are transmitted to Fleet Numerical Oceanography Center 
(FNOC) by approximately 3 h after 00/12 GMT. These 
soundings are used tc generate the 12-hourly forecasts of 
the 500 mb heights (in this case) . These generated fields 
would be available tc Joint Typhoon Warning Center (JTWC) 
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Guam prior to the next warning at approximately 6 h 
after 00/12 GMT. The computer program to decorapos-e the 
D-value fields (deviations from the standard atomsphere 
gs opoten tial) into the eigenvectors (EOF coefficients) and 
compute the persistence statistics, would take at most a few 
minutes, and could actually be done on a hand held 
calculator (Shaffer, 1982). This would make the warning 

time increments equal to 30, 4 8 and 66 h if the forecast, 

technique makes use of the analyzed and forecast 
gsopotential height fields. This time delay problem could 
be partially eliminated by using the previous 12-h forecast 
fields to define steering, which would permit the technique 
to be available at the time JTWC prepares their forecast at 
00/12 GMT. It is believed that little effect would be 
apparent for the shorter forecast periods, although there 
would likely be a larger impact on the extended time period 
forecasts. 
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II. EOF ANALYSIS OF SYNOPTIC FORCING APPROACH 
TO~2ECIJR71TURE 



This research is a sequel to a synoptic map type scheme 
of Brown (1981) and an EOF representation of the synoptic 
forcing to predict tropical cyclone motion in the western 
North Pacific Ocean by Shaffer (1932) . The data fields and 
selection of storm criteria (see descriptions by 3rown, 
1981; Shaffer, 1932; and Shaffer and Elsberry, 1932) are 
exactly the same in all three studies. A short summary of 
the fields and criteria follow. 
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Figure 1 . 



Grid Arrangement [Jsei in Statistical 
-Climatological Track Prediction. 
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The 500 mb D-values are extracted on a movable 120-point 
grid (Fig. 1) over 35° in latitude and 70° in longitude with 
a 5° increment. The tropical cyclone is always centered at 
grid point 70 which is 10° from the southernmost boundary 
and 40° from the western boundary. This choice is made :o 
include both a portion of the aid-latitude flow 

characteristics and a smaller section of the more uniform 
tropical height fields to the south of the storm. The 
movable feature is incorporated to minimize the seasonal 
variation that would occur if the grid was stationary. The 
earth’s curvature (convergence of the meridians) is not 
included over the grid but could be easily added for 
operational use. The D-values on the standard 63 x 63 grid 
from the FNOC archives are deviations from the standard 
atmosphere geopotential and are fitted to the 120 point grid 
by a Bessel linear interpolation scheme. 

Requirements for tropical cyclones to be included in 
this study are: 

* The tropical cyclone must b e of at least trooical storm 
intensity (>17 m/sec) and be located to the west of 
180°; 

• The cyclone must exist at tropical storm intensity for 
at least 30 hours as analyzed by JTWC (Guam) ; 



• The cyclone must be located between 10°M and 25°N' (This 
is done to prevent the grid from entering the southern 
hemisphere and to avoid oossible problems arising from 
the convergence of the meridians at higher latitudes) ; 
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• The storms must exist at 00/12 GMT for the scheme to be 
able to use the corresponding FNOC data analyses; 

• Each storm observation included must have 35 hours 
between a prior observation of the same storm to insure 
pseudo-independence within the sample. This very 
important restriction insures that the statistical 
results are unbiased. 

One approach to assess spatial variations in geophysical 
fields is through the use of Empirical Orthogonal Function 
(EOF) analysis (also known as Principal Component analysis). 
EOF analysis may be applied to any multi-dimensional data 
set to separate signal from noise, and identify significant 
spatial correlation patterns (Hotelling (1933); Kutzbach 
(1967); Morrison (1967); 3arnatt ( 1978); Richmac ( 1980); 
Shaffer (1982) and others). For geophysical use, the data 
are placed on the moveable grid with each grid point 
representing a dimension in space. The data are represented 
as an MxN matrix, with M the number cf grid points, and N 
the number of data cases. Each column of this matrix 

represents a single realization of the observed geophysical 
field. An MxM matrix R is then formed by multiplying the 
data matrix by its transpose. The form of R depends upon 
the form in the initial data matrix, but is generally, 
either a covariance or correlation matrix of the grid 
points. If one wishes to give equal weighting to each grid 
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on the other 



point, the correlation matrix is formed. 1 If, 
hand, the desire is to give greatest weighting to those 
dimensions (grid points) with the largest univariate 
variation, the covariance matrix is formed. 

The choice of correlation versus covariance matrix is 
critical for meteorological fields which cover a large 
spatial domain. For example, in analyzing geopotential 
fields, mid-latitude grid points have much greater 
variability than tropical grid points. If the aim is to 
focus on mid-latitude features, the covariance matrix form 
of R is preferable. If, on the other hand, the wish is tc 
determine variability patterns over the entire domain, the 
correlation matrix form of R is favored. The correlation 
matrix approach mutes the EOF response in higher variability 
regions (mid-latitudes), while amplifying the response in 
regions cf lower variability. Both approaches were tested 
in this study. With the covariance matrix approach, the 
total explained variation in the D-value fields for the 
first 10 eigenvectors was more than 9U&, which was 9 % 
greater than the correlation eigenvectors. Regression 
analyses from the covariance method were shown to be 
inferior to these from the correlation method. Thus, the 
correlation matrix method is used in this study. 
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After forming the matrix a, minimization of the total 
variability leads directly to the eigenvalue problem 
(Morrison, 1967). The eigenvector associated with the 
largest eigenvalue is the single vector, out of all possible 
vectors in space, that minimizes the summed variation 
(squared) of all of the observation vectors in M space. 
Furthermore, this eigenmode explains exactly 



x, 



M 

E X, 



i-1 



of the total variation in the observations, where 



e X is t 



he 



eigenval ue 


of the ith mode. 


Since the 


trace of 


correlation 


matrix is identically 


equal to the 


order of 



matrix, the percent of total variation explained by any 
eigenvector is simply 

X, 



s X, 



The elements of each eigenvector may be loosely interpreted 
as the multiple correlation of the corresponding grid point 
with all other grid points, simultaneously. The exact 
meaning of the eigenvector elements (loadings) is given by 
Morrison ( 1967) . 
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The aim of this project is to investigate spatial 
variation structures of analyzed 500 mb D-value fields for 
the western North Pacific Ocean. Therefore, we choose the 
correlation matrix approach. The initial step in data 
processing is to point-nor malize the data. Specifically, we 
use 



Dbar (i) 



_L E Dti,j) 

N j-i 



D*<i, j) = D (i,j) - Dbar(i) 

N 



[£■- 

l -- 1 



D (i, j) - Dbar (i) ) **2 2 



N-1 



P 



where j is the map (case) index and i is the grid point 
index. Thus, D* (i,j) is simply a normally distributed 

random variable, with a mean at each grid point over all 
maps of zero and a variance of 1. This allows an 
equally- weighted, grid point interpretation of rhe 

eigenvectors. For this study, rhe daoa analysis was 
completed using the 500 mb fields only. 

The data consist of D-values of geopotential, (the 
difference between the observed geopotential and the 

standard gecpotent ial, so that a positive D-value 
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corresponds to higher than standard geopotentials ) , m 
meters. EOF analysis was performed for 120-dimension space 
(8x15 grid) in Fig. 1. Examples of the first 9 EOF 

eigenvectors (correlation case) are shown in Fig. 2. It 
must be stressed that the negative of the pattern shown in 
these figures is also a solution to the problem. (For 

examples of 700mb and 850mb eigenvectors see Shaffer 
(1 982) .) 




Figure 2. Eiaenvectors 1-9 ;iuj.t lplied, by 100, The X, 

Indicares the Location or the Tropical. Cycior.e, 
Grid Point 70. 

The number of eigenmodes to retain in each case is 
determined by comparing the eigenvalues of the analyzed 
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geo potential fields with the eigenvalues of a Monts Carlo 
simulation of random fields as used by Preisendorf er and 
Barnett (1977). The eigenmodss retained are those that are 
larger than the random eigenvalues plus 3 standard 

deviations. By using this approach we are 99% confident 

that the only eignvectors retained are those containing 
actual signal, and that we are discarding noise only. Using 
this selection procedure, the actual number of significant 
eigenvectors (SOF's) retained is limited to 10. This number 
of EOFs explains approximately 85% of the variation in the 
fields. Table I outlines the increase of explained variance 
as the number of significant eigenvectors retained is 

increased. Fig. 3 is the graphical representation of Table 

I. 

The total sample of the 504 storms during the 1967 to 
1976 time period is divided into a dependent sample of 454 

cases and a 50-case independent sample for use as the test 

sample. This is the identical sample of 50 randomly 

selected independent storms that Shaffer (1982) tested for 
displacement forecasts. Shaffer (1982) and Shaffer and 
Eisberry (1982) showed that for this independent sample of 
tropical cyclones, the forecast error from a technique based 
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TABLE I 



Eigenvalues and Cumulative Pe 
of the Normalized D-value Fields 


rcent Explained Variance 
at 500 mb. (Shaffer# 198 


Mode 


Eigenvalue 


Explained Variance 


1 


39.36 


33 


2 


21 .60 


5 1 


3 


9 .29 


59 


4 


7 .43 


65 


5 


6 .08 


70 


6 


5 .29 


75 


7 


4 .00 


78 


8 


3.13 


8 1 


9 


2.49 


83 


10 


2.11 


85 


15 


1 .03 


91 


20 


.53 


94 


4 0 


.11 


98 


60 


.04 


99 


120 


.00 


100 


on an 


EOF represen tar ion of 


zhe syaopzic forcing 



decreased by an average of 17% * hen compared z o -he official 
JTWC (Guam) forecast error. 
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Figure 3. Eigenvalues of the 500 mb D-value Fields 
(ctrcles) and the don te Carlo generated 
Eigenvalues (triangles) Based on 100 Simulations 
of Random Data. 
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• DEFIN ITIO N OF RECUR VAT UR E 

The definition of a recurvature case can be thought of 
as a set of restrictions based on the observed track of the 
tropical cyclone. The definition of recurvature that is 
used ir. this research is than described by Leftvich (1978, 
1982): 

"Recurvature is a net displacement northward of 315° 
during the forecast period or the attaining of 

northeastward morion during the 12 hours prior to the end 
of the forecast period". 

Storms are not considered if their past 12 hour ■'motion is to 
the east. A graphical descriprion of recurvature and 
non-recu rvarure cases is given in Fig. 4. 

Some of the cases present rarher difficult forecast 
situations. Figure 5 depicrs Typhoon Trix (initial position 
at OOGMT 21 August 1971) , which for the initial 36-h period 
would be considered to have recurved, whereas for the 54- 
ana 72-h period it would be considered as a non- recurvar ure 
case. Figure 6 depicts Typhoon Polly (initial position at 
12GMT 29 August, 1974) , which for the 36-h period is 
considered to have recurved, but then moves westerly to a 
54-h position which is considered as having not recurved. 
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and finally moves north-north* est to be considered as a 
recurvature case at 72 h. Typhoon Cora (initial position of 
12GMT 02 October 1975) shown in Fig. 7 is an example of a 
recurvature case at all three periods, whereas Typhoon 
Alice (initial position at 12GNT 16 September 1975) in Fig. 
8, is a non-recurvature case for all three periods. Typhoon 
Anita (shown in Fig. 9) would net be considered for the 
12GMT 22 July to 00GM1 23 July initial positions due to the 
eastward motion of the tropical cyclone, but would be 
considered for the 12GMT 23 July to 00G11T 24 July period. 

A summary of the recurvature statistics for the 504 
storm sample is given in Table II. The 6 h recurvature 
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Figure 5 



Figure 5 




. Typhoon Trix Case Illustrating an Apoaren 
Recurvatur-s Case a* 3 6 h, and Nor.-recurva 
54 and 72 h 




25 



C» fl 



Fi gure 



Figure £ 

st atisti 
ea rlier. 




. Typhoon Cora Case Illustrating a Common 
Recurvature Case 




cs were obtained in the same manner as def 
except, that the motion during the past 6 h was 



t u r e 

inea 

used 
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Figure 9. Typhoon Anita Case Illustrating the Exclusion o 
the 12GMT 22 July to OOGMT 23 Julv Periods Due 
the Eastward Motion 

for the second consideration for recurvature vice the past 
12 h motion. All past meridional and zonal motions were 
obtained from the warning positions of each individual 
storm. This was done to simulate actual operational 

conditions (Shaffer, 1982). 

A graph of the recurvature statistics shown in Table II 
is presented as Fig. 10. There is a slow increase of 
probability of recurvature with time. This increase may be 
explained by the following reasoning. A tropical cyclone in 
its developmental stages is more often in the easterly trade 
wind regime and moving west of northwest (non-recurver as 
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1 1 r-t> 



TABLE II 



Number of Recurving (rev). Non-recurving (r.rcv) , N<? History 
(nohist) , Not Considered (notcon) and the ProbaDility of 
Recurvature (prob) for the 504 Tropical Cyclone Sample 



tau 


rev 


nrev 


nohist notcon prob 


06 


131 


234 


0 


139 .359 


12 


125 


240 


0 


139 .342 


18 


130 


235 


0 


139 .356 


24 


134 


231 


0 


139 .367 


30 


135 


230 


0 


139 .370 


36 


158 


196 


22 


128 .446 


42 


1 19 


169 


1 15 


101 .413 


48 


125 


163 


1 15 


101 .434 


54 


133 


154 


1 16 


101 .463 


60 


133 


139 


1 42 


90 .489 


66 


96 


1 11 


2 30 


67 .464 


72 


103 


104 


2 30 


67 .498 


78 


107 


100 


2 30 


67 .517 


34 


99 


87 


2 54 


64 .532 


defined) . As the tropical cy 


cl on> 


e intensifies and moves to 


higher latitude 


S/ 


it beg 


in s 


to interact wirh the 


mid-latitude wes 


ter lies and 


th er< 


efore is more likely ro 


recurve. Cheng- 


Lan 


and Sadler 


(1982) describe the racsr 


frequent track of 


a r 


ecur ving 


tropical cyclone as "a regular 


clockwise parabola". 


That is. 


a s ■ 


the tropical cyclone moves 


westward it also 


move 


s norzhwa 


ri. 


wirh a greater likelihood 



of being north of a 315° line from the initial position. 
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6.0 12.0 18.0 24.0 30.0 36.0 42.0 48.0 54.0 60.0 66.0 72.0 73.0 64.0 

TIME IN HOURS 

Figure 10. Probability of Recur vature (504 Case Sample 
Climatology) 

The statistics contained in Table II and Fig. 10 are 
designated as the sample climatological probabilities of 
recurvature. If should be mentioned that the independent 
sample statistics for recurvature for all time periods are 
within 7% of the sample average probability, which indicates 
consistency between the two samples. These values will be 
used later as a standard for evaluation of various forecast 
schemes in predicting recurvature at specified time 
intervals. That is, a scheme will be said to have skiII if 
it predicts recurvature more accurately than the 
climatological probability. 
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IV. REGRESSION ANALYSIS 



Neumann and Lawrence (1975) demonstrated that a 
statistical-dynamical model using prognostic and 
persists nce-related variables in a nonlinear regression 
scheme shewed much promise in forecasting Atlantic ocean 
tropical cyclone motion. Shaffer (1932) also used a 
st atisri cal- dynamica 1 scheme that used 10 30F coefficients 

to represent the geopetential field (analysis time only) and 
a linear regression scheme to predict tropical cyclone 
motion. The 30F analysis not only simplifies the 
representation of the synoptic forcing, it also reduces 
significantly the cos: of producing an accurate and viable 

f orecast . 

The importance of the nonlinear rerms in regression 
schemes has been demonstrated by Neumann and Lawrence (1975) 
and Paegle and Haslam (1982). Although nonlinear terms were 
included in those studies, economic considerations and the 
simplicity of application of a linear model led to ins 
adoption in this study. 
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To test the effect of including ED F coefficients as 
predictors of tropical cyclone motion, three separate trials 
for each of the prescribed times were completed. The first 
test (refered to as PERS) was accomplished using only the 
persistence-related parameters and the initial latitude and 
longitude. A description of all potential predictors is 
contained in Table III. (All latitude and longitude 

displacements (displ) are in nautical miles.) The second 
experiment (referred to as EOFS) utilized 10 EOF 
coefficients and the initial latitude and longitude. The 
third experiment used a combination of the persistence and 
EOF predictors and is noted as ?&S. It should be mentioned 
that the cyclone maximum sustained wind speeds at 12, 24 and 
36 h prior to the initial observation were also incorporated 
as possible predictors. This was dona in accordance with 
Riehl's hypothesis (1971) of maximum intensity being 

associated with recurvature for Atlantic tropical cyclones. 
However, none of these additional predictors were chosen 
during the stepwise regression. 

The 504 case tropical cyclone sample was analyzed for 
recurvature as described in Chapter III. A value of 1.0 was 
assigned to the storms that recurved and a value of 0.0 was 
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TABLE I II 



Potential Predictors Used to 
Equations for Each of the 


Develop the Regression 
Three Time Intervals 


Potential Predictor 
Variable Number 


Name 


Descrj.pt ion 
(Coef ici ent-Coef ) 


1 


CO FI 


Eigenvector 1 

Coefficient 


2 


COF2 


Eigenvector 2 

Coefficient 


3 


COF3 


Eigenvector 3 

Coe f f icient 


4 


COF4 


Eigenvector 4 

Coefficient 


5 


C0F5 


Eigenvector 5 

Coefficient 


6 


C0F6 


Eigenvector 6 

Coe f f icient 


7 


COF7 


Eigenvector 7 

Coefficient 


8 


C0F8 


Eigenvector 8 

Coefficient 


9 


C0F9 


Eigenvector 9 

Coe ff icient 


10 


CO F 1 0 


Eigenvector 10 

Coef ficient 


11' 


PLATO 


Initial Cyclone 
Latitude 


12 


PL AT 1 


Past 12 hour 
meridional displ 


13 


PLAT2 


Past 24 hour 
meridional displ 


14 


PLAT3 


Past 36 hour 
meridional disol 


15 


PLAT4 


Meridional displ 
for prior 12-24 hrs 


16 


PLAT5 


Meridional displ 
for prior 12-36 hrs 


17 


PLAT6 


Meridional displ 
for prior 24-36 hrs 


18 


PLONO 


Initial cyclone 
longitude 


19 


PLON1 


Past 12 hour 
zonal disDl 


20 


PLON2 


Paso 24 hour 
zonal displ 


21 


PLON3 


Past 36 hour 
zonal displ 


22 


PLON4 


Zonal disci 
for prior 12-24 hrs 


23 


PLON5 


Zonal displ 
for prior 12-36 hrs 


24 


PLON6 


Zonal displ 
for prior 24-36 hrs 


assigned to the cycl 


ones rhat had not recurved at the time 


of interest. A value 


of 999. was assigned to cyclones that 


had no history at the 


time period 


. of interest. 
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The BMDP stepwise linear regression program screened the 
dependent data set for storms with past 12 h longitude 
displacement (PL0N1) that was to the east and also tropical 
storms that had no history. This screening accounts for the 
different number of cases considered for the regression 
routine and the cliiratology cases at the specific time 
periods of interest (see Table II). The equations at 36, 54 
and 72 h (see Tables 17, V, 7 1 below) are derived from a 

stepwise linear regression using an F - ratio of 4.0. The 
F-ratio is a test of the significance of the independent 
variable coefficients in the regression equations. The 
number of cases of the dependent sample that are selected 
for each regression period is noted in the corresponding 
table. The number of cases decreases with time due tc the 
restrictions for inclusion discussed in Chapter II. Of note 
is that most of the cases at 72 h are a subset of the 54 h 
data set, which is in turn a subset of the 36 h data set. 
The R Squared Statistic ,RSQ, (percent of variance in the 
predictar.c explained by the regression equation) will be 
used for discussion of the results. 



33 



The linear regression equations produced from the 
potential predictors of recurvature or non-recurvatur e are 
of the form: 



? (0: 1) = aO + a 1X1 + a2X2 + a3X3 + 



+ a7X7 



Tables IV, V, and VI list the coefficients and variables 
for each of the three tests generated from the dependent 
data sets. The number of cases used to generate the 
coefficients are 259, 206 and 133 respectively. 

TABLE IV 

Regression liquations for 36 h Recurvaturs for Each of the 
Three Methods Tested (Dependent Sample) , and Mean Values 
and the Standard Deviation from the Mean (STD) 



PERS EOFS 


PSE 


MEAN 


STD 


intercept . 1852 -^.2871 


.5502 






PLATO .0167 




19.4 


5.2(0) 


PLONO .0054 




133.2 


12.9 (°) 


PL A T 2 .0017 


. 00 14 


96.4 


90. 2 (NM) 


PLON1 .0031 


. 0031 


-92.6 


55.3 (NM) 


COE 1 -.0326 


-.02 15 


.81 


5.12 


COF 3 .0367 




-.24 


2.89 


COF 5 -.0302 




. 20 


2.34 


CO F 7 .0341 


. 02 99 


-.07 


1.91 


COF 8 .0493 




-.20 


1.77 


COF 10 -.0545 


-.03 97 


-. 10 


1.48 


Each regression coefficient 


is unique 


except for 


the valu 


of PLON1 fcr the PERS and PSE 


test at 


36 hours. 


It shoul 
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TABLE V 



Regression 


Equations 


for 54 h Recurvatur 
Table IV 


e as Descri 


bed in 




PERS 


EOFS 


P5E 


MEAN 


STD 


INTERCEPT . 


4446 


.4957 


.4873 






PLAT2 


0019 




. 00 14 


95.6 


86.6 


PL0N1 


0022 




.00 16 


-93.0 


56.6 


COF 1 




-.0312 


- -. 0252 


1.23 


4.3 


COF 3 




.0371 


.0237 


-. 53 


2.74 


COF 8 




.0618 


. 0465 


-.21 


1 . 62 


COF 10 




-.0514 




.05 


1.50 






TABLE VI 






Regression 


Equations 


for 72 h Recurvatur 
Table IV 


e as Descri 


ned in 




PERS 


EOFS 


PSS 


MEAN 


SID 


intercept 


5487 


.5311 


. 5821 






PLAT4 


0028 




. 0026 


44. 1 


48.1 


PLON1 


0027 




.00 24 


-89. 1 


52. 2 


COF 1 




-.0 3 24 


-. 0230 


2. 02 


4.55 


COF 3 




.0410 




-. 75 


2.77 


be noticed that the value of 


PL ON 1 must 


n ecessarily 


be less 


than zero. 


because 


a PL ONI 


greater than zero represents 


easterly mov 


ement , 


which is 


c cr.sidered 


as having 


already 


recurved . 












The total number of combinations 


of predictors and 


regression coefficients will 


not be discussed, but 


several 


combinations 


will be 


mentioned. PLAT2 


(past 24 h 


latitude 
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displacement) is almost always greater than zero (the 
exception being trochoiaal oscillating storms), and makes a 
considerable contribution towards recurvature (predictand 
close to 1.0) especially for large displacements. PLATO and 
PLONO (initial latitude and longitude) are always greater 
than zero, with larger values being favorable for 
recurvature prediction. The coefficient of the first 
eigenvector (COF 1) has the opposite effect. Physically, 
the stronger the zonal pressure gradient (large COF 1), the 
lower the probability of recurvature. The opposite is also 
true; the probability of recurvature is higher for 
extremely weak zonal pressure gradients as illustrated in 
Fig. 11. 

Analysis of the sample mean values of the coefficients 
of EOF 1 (assuming a normal distribution about the mean) 
indicates that the number of recurving tropical cyclones 
would be less than the number of non-recurving storms. For 
example, the number cf recurving storms at 36 h (Table II) 
is 158 and the number of non -r e curv ing storms is 196. The 
coefficient of EOF 5 is also quite interesting. The more 
the EOF 5 pattern describes the synoptic forcing (the larger 
the value of COF 5) , the more likely the storm will not 
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Figure 1 1 
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Displacement (Shaft ar, 1982) 
recurve. The converse is also true in that the stronger the 
inverse or the negative of the pattern (a strong shortwave 



between two ridges) , the more likely the cyclone is to 
recurve. 



This 
possible 
and their 



is by no 
54 linear 
inverses 



means a complete discussion of all the 
ombinations of predictors, coefficients 
It does show how a simple combination 
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of EOF patterns can be applied with straight-forward 
meteorological reasoning. 

A. 36 H REGRESSION CASE ANALYSIS 

As expected, the predictors selected for the PE3S 
regression analysis included the past 12-h longitude 
displacement, the past 24-h latitude displacement and the 
initial latitude. Simple extrapolation techniques do well 
for these relatively short-time forecasts. The longitude 
displacement over the past 12 hours reduces the variance in 
the predictand (increases RSQ) slightly more than the 
latitude displacement. The selection of initial latitude as 
a predictor is consistent with the variation of the 
probability of recurvature with time (and therefore latitude 
for an idealized storm track) as discussed in Chapter II. 
The total explained variance of the predictand for the 
persistence-related predictors was 29.2^. 

The individual EOF predictors reduce less of the 
variance than the persistence-type predictors (Table VII). 
Notice that eigenvector 1 (which is an indicator of the 
strength of the zonal pressure gradient and large global 
wave number) causes the largest reduction in the variance of 
the predictand. The remaining five eigenvectors that were 
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selected are indicative of smaller scale forcing features, 
and reduce the variance of the predictand much less that the 
first eigenvector. Shaffer (1982) indicated the possibility 
of weak northwest motion (or southeast motion for the 
inverse mode) would he associated with eigenvector 7 and a 
strong meriodional forcing being associated with eigenvector 
3. Of special interest is eigenvector 5. Although the 
reduction in variance associated with this eigenvector is 
small, its inverse pattern appears similar to a strong short 
wave in the westerlies which could advect the storm 
northward. The last predictor selected, the initial 
longitude, again is consistent with the variation of the 
recurvature probability with time. The total of the 

explained variance for the 30?S test was 26 . 5 % , which is 

2.7% less than using persistence only predictors. 

The selection of predictors for the final test shows a 
combination of both PERS and S3 FS predictors. The largest 
reduction of the variance occurs with past motion 

(persistence related) parameters. The strength of the zonal 
pressure gradient appears as the next most important 
predictor with smaller scale features (EOFS 10 and 7) 
contributing small increases to the explained variance. The 
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total of the explained variance of the predictand is 34.0% 
which is greater than either the PERS or EOFS methods, but 
still only a small fraction of rhe total variation. 

Part of the reason than more of the predidictand 
variation could not be explained may be related to the 
arbitrary definition of recurvature, e.g. a movement no 

320° is recurvature (value 1.0) whereas 310° is not (value 
of 0.0). When a yes/no decision is applied across a line as 
in this situation, one cannot expect an explained variance 
near 100%. The method could still provide useful guidance 
for the forecaster if the percentage of correct forecasts is 
higher than other methods available. 

The entry of PLATO, and than the subsequent delerion as 
the EOF information is added, is an in dican ion of 
overlapping information. Notice (Table '/II) that PLATO does 
not enter in scheme (c) at all. 
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TABLE 7 1 



Summary of Stepwise Linear Regression Results for 36 h (259 
cases) R Squared (RSQ) , Increase rn R Squared (RSQ Incr) 
and F to Enter or Remove 

a) PEESistence only 

F remove 



Step 




Parameter 


RSQ 




RSQ 


Incr 


F 


9 


nter 


1 




PLON1 


151 




. 1 


51 




45 


.65 


2 




PLAT2 


263 




. 1 


12 




39 


.06 


3 




PLATO 


292 




.0 


29 




10 


.42 


b) 


EOFS 


only 
















1 




CGF1 


121 




. 1 


21 




35 


.40 


2 




PLATO 


161 




.0 


40 




12 


.04 


3 




COF 8 


187 




.0 


26 




8 


.20 


4 




COF7 


206 




.0 


19 




6 


.04 


5 




COF3 


227 




.0 


21 




6 


.85 


6 




COF10 


24 3 




.0 


16 




5 


.47 


7 




PLATO 


234 




-.0 


09 








8 




COF5 


248 




.0 


14 




4 


.54 


9 




PLONO 


265 




.0 


17 




5 


.75 


c) 


Pers 


istence & EOF 


S 














1 




PLON1 


151 




.1 


51 




45 


.65 


2 




PLAT 2 


263 




. 1 


12 




39 


.06 


3 




COF1 


313 




. 0 


50 




18 


.54 


4 




CCF10 


327 




.0 


14 




5 


. 10 


5 




COF7 


340 




. 0 


13 




5 


.18 


5. 


54 H 


REGRESSION C 


ASZ 


ANALYST 


S 










The 


selection of 


PERS 


pr; 


rdic 


tors at 


54 


b 


(Ta 


the 


same 


as at 36 h 


(Table ’ 


/I) 


except 


f 0 


r 


the 


the 


init 


ial latitude. 


H 


ad i 


: he 


F-test 


bee 


n 


set 


4.0, 


t 


he initial 


lat 


itude 


wo ul d 


hav 


o 


be 


subseque 


nt to PL ON 1 , 


as 


was 


the 


case f 


or 


36 


■h. 


lat i 


tudi 


nal dispiacem 


ent 


for th 


e past 


24 




h de 


var i 


ance 


of the pre 


diet 


and 


mor 


e than 


a 


i d 


the 


iongitud 


e predictor. 


The 


se i 


: wo 


predict 


ors 


c 


ompr; 



2.94 



h Cl 
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the variance of the predictand, which is almost half the 
value obtained at 36 h. This decrease is not surprising, 
because simple extrapolation techniques have less validity 
at 54 h than at 36 h. 

The number of predictors selected in the EOFS test 
decreased to four, compared to seven at 36 h. These four 

predictors at 54 h are all included within the first five 
predictors at 36 h. The north/south pressure gradient (EOF 
1) again contributed to the largest reduction in the 
variance of the predicta nd (see Figure 2). EOF 8 appears to 
be similar to EOF 5 in that it is associated with a 



mid-latitude short 


wave comb 


ined with 


an area 


o f 


h 


igh 


pressure directly 


south of the 


tropical 


storm . 


EOF 


6 


is 


correlated directly 


with recur 


vature as 


are EOFS 


1 


and 


3. 



EOF 10 is indirectly correlated with recurvature due to the 
high pressure area to the northwest of the storm center. 
The stronger this pattern, the less likely it is that the 
storm will recurve. The total of the explained variance of 
the predictand is 19.4^ which is smaller than the value for 
36 h using EOFS only, but it is larger than 54 h regression 
using persistence parameters only. This result suggests 
that synoptic parameters begin to explain more of the motion 
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at later times than does the extrapolation associated with 



persistence variables. 

The combination of the two types of predictors explain 
24.65? of the variance of the predictand. The only EOF 
predictor that does not appear in P&E is eigenvector 10. 
The inclusion of the eigenvectors as predictors increased 
the RSQ amount by 8.2ft. The decrease from 34.05? explained 

variance at 36 h to 24.6ft at 54 h was expected due to the 
likelihood of more tropical cyclone interactions with the 
synoptic environment. 

TABLE VIII 

Summary of Steowise Linear Regression Results for 54 h 
(206 cases) as Described in Table VII. 



a) FERSistence only 



Ste p 


Parameter 




RSQ 


RSQ Incr 


F ei 


iter 


1 


PLAT2 


• 


102 


. 1 02 


23. 


. 26 


2 


PLON1 


• 


164 


.0 62 


15. 


, 06 


b) 


ZOFS only 












1 


CO FI 




C74 


.0 74 


16 . 


.20 


2 


COF8 


• 


127 


.0 53 


12. 


. 42 


3 


COF3 


• 


170 


.043 


10, 


.41 


4 


COF10 


• 


194 


.024 


6. 


. 03 


c) 


Persistence & 


EOF 


S 








1 


PL AT 2 


• 


102 


. 1 02 


23. 


. 26 


2 


PLON1 


• 


164 


.0 62 


15. 


.06 


3 


CO FI 


• 


203 


.039 


9. 


, 93 


4 


C0F8 


• 


230 


. 0 27 


6 . 


.97 


5 


COF3 


• 


246 


.0 16 


4. 


, 15 
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C. 72 H REGRESSION CASE ANALYSIS 

The selection of predictors for the 72 h prediction 
includes a new predictor, PLAT 4 (latitude displacement fcr 
the period 12 to 24 h prior to the initial time). This 
predictor is inherently positive and therefore adds to the 
probability of recurvature. The past 12 h longitude 
displacement is the remaining predictor chosen. The total 
of the explained variance using these two predictors is 
15.0??. Although this is less than half of the value of the 
36 h regression, it is only slightly less than the 16.5% 
obtained for the 54 h regression. 

TABLE IX 



Summary of Stepwise Linear Regression Results for 72 h 
(138 cases) as Described in Table VII. 



a) 


PERS 


istence 


only 














Step 




Paramet 


er 




FSQ 


RSQ 


Incr 


F er 


iter F 


remove 


1 




PLGN1 




• 


074 


.0 


74 


10. 


, 84 




2 




PL AT 4 




• 


150 


.0 


76 


12. 


,01 




b) 


EOFS 


only 


















1 




CO FI 




• 


076 


.0 


76 


1 1 . 


, 27 




2 




C0F3 




• 


123 


. 0 


52 


8. 


.03 




c) 


Pers 


istence 


& 


EOF 


C 












1 




COF1 




• 


077 


. 0 


77 


1 1 . 


.27 




2 




PLON1 




• 


132 


. 0 


55 


8. 


.63 




3 




PLAT4 




• 


152 


. 0 


60 


9. 


, 94 






The 


EOFS-cnly 


re 


gression 


or 


oduced re 


sul* 


:s that 


explain 


less 


of 


the 


var 


ian 


ce of 


the 


predict 


a nd 


than d 


id the 
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persistence-type predictors. Coefficients 1 and 3 explain 
12.8? of the total variance of the predictand, which is 
almost identical to the value for the first steps of the 54 
h regression. However, two different coefficients entered 
the 54 h regression equation. 

The combination cf two groups of predictors is again 



superior to 


either individual 


group. At 72 h 


the 


north/so ut h 


pressure gradient 


is 


more important than 


the 


longitudinal 


displacements 84 


h 


earlier (PL0N1) or 


the 


latitude dis 


placement 84 to 


9 6 


h prior to the ini 


tial 



position (PLAT4) . Use of extrapolation-type methods can not 
be expected to contribute as much as the synoptic parameters 
at such long time intervals. As the tropical cyclone moves 
into the mid-latitudes, persistence is comparatively less 
useful as a predictor than in the tropics. 

The amount of explained variance of the racurvature 
event decreased with time for each of the three test groups. 
The superior method (as far as explaining the variability of 
the predictand) is the combination of persistence-related 
variables and the EOF coefficients which describe the 
synoptic patterns. Although the PS E technique consistently 
explained more of the variance than the other two methods. 
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further zesting is necessary 
predictions are superior to 
climatology or mere chance. 



to determine if 
other techniques 



the 


EOF 


such 


as 
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V 



. TESTING THE DEGRESSION ANALYSIS 

Choosing a forecast verification scheme is very 
important for testing the ability of a regression model or a 
prognostic method tc forecast a specific event. The 
geophysics community is just one of the many scientific 
communities that desire an accurate and acceptable method 
for describing the skill of a given model. Muller (1944) 

provided a literature survey of 55 short-range weather 
verification methods that were developed during 1884 to 1943 
in six different countries. Muller divided the schemes ir.ro 
two areas; comparisons of the forecasts and observations as 
a ratio (percentage correct) and verification methods which 
are compared to chance or climatology, the so-called 'blind' 
forecast s. 

In the early verification studies, the schemes for 
comparing forecasts with observations had to be simple for 
the ease of application. Verification methods today compare 
the prediction skill against some climatological value or to 
pure chance as a point of reference. 
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Three skill scores will be discussed: the Heidke Skill 

Score (HSS), the Brier P-Score (BPS), and the percent of 
correct forecasts (PCF) . The Heidke skill score can be 
written as: 

HSS = VP -C 
3? - C 

where V? is the number of correct or verifying forecasts of 
an event (either correctly forecast recur vers (CFR) , i.e., 

predictand is greater than or equal to an assigned 
climatology probability (Table 1 1) ; or a correctly forecast 
non-recurver (CFNR) , i.e., the predictand is less than the 

climatology probability) , T is the total number of forecasts 
made and C is equal tc the number of correct forecasts that 
would be expected from chance. In this test C is always 
one-half of the total number of forecasts. The HSS would 
have a maximum value of 1.0 for perfect forecasts and a 
value of zero when the number of correct forecasts is equal 
to the number of forecasts expected to be correct by using 
chance. Negative values occur if the method is inferior to 
the comparative methods. 
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The Brier or P-Score has been widely accepted as a 
method to test probability forecasts. The binary form of 
the EPS can be written as: 



where FCST is the regression equation predictand and ACT is 
either 1.0 if the storm actually recurved or zero if 
recurvature did not occur. A perfect BPS is a value of zero 
and the worst possible forecast wculd be indicated as 2.0. 
It is readily apparent that the BPS is actually the mean 
square error of the set of probabilistic forecasts. 

The percent of correctly forecast events (PCF) can be 
written as: 



where CFR, CFRR and T are the same variables as described 
earlier. 

As discussed in Chapter IV, the predictand values should 
lie cn the continuous interval between 1.0 and 0.0. The 
unfiltered predictand statistics for both the dependent and 
independent samples for all periods appear in Table XI. 



N 




ACT (i) ) **2 



PCF = CFR + CFNR 



T 
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Values larger than 1.0 (indicative of a strong recurvature 
signal) and smaller than 0.0 (indicating strong 

non-^recu rvature signal) exist as maximum and minimum values. 
These values outside of the defined range are due to the 
assumed linear relationship between the predictors. Table 
XII provides similar statistics for the predictahds after 
values that were more than 1.0 and less than 2.0 were set 
equal to 1.0 and negative values were sat equal to 0.0. 
Comparison of the dependent and independent samples of both 
tables indicate that the two sets are quite similar. The 
different number of cases included in several tests is due 
to the predictors PLAT2 and PL0N1 being excessively large# 
(purposely assigned flagged values) which causes several 
predictand values for PERS and P&E to be in excess of 2.0, 
the upper filter limit. 

Tables XI and XII indicate that the variance of the 
dependent predictand is greater than or equal to the 
variance of the independent predictands in 15 of 18 
comparisons. This is evidently due to the small size of the 
independent sample. 

It also should be noticed that the dependent sample 
variances are largest at 36 h and decrease in time. In all 
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dependent data set comparisons, the variance of the P&E tesr 
is larger than the other two tests. This is most likely due 
to the 'signal* of either recurvature or non-recurvature 
being strongest for this set which includes both persistence 
and synoptic features. The smaller values of the variance 
in Table XII, as compared to Table XI, is due to the 
filtering of the predictands, which is described above. 

TABLE XI 

Predictand Statistics for the 36, 54 and 72 h Before 
Filtering. The Number of Cases, Mean, Maximum and Minimum 
Values (max.min) and the Variance (var) of the Predicted 
Values ror the Dependent and Independent Data Sets are Shown 
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TABLE XII 
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Table XIII mst be carefully interpreted. when using “he 
HSS, any predictand that is greater than or equal to a 



52 



specific value (in this case a climatological percentage of 
recurvature for that forecast interval) is given a forecast 
recurvature value of 1.0. Any predictand that is less than 
the climatological value is considered to be a 

non-recurvature case, and is assigned a value of 0.0. These 
'newly assigned' predictands are then compared to the actual 
storm recurvature values. If the forecast ana the actual 
value match, a verifying forecast is tallied (CFS or CFNR is 
incremented by one). 

A graphical comparison of HSS versus a varying 

climatological value cf probaoility of recurvatura (used for 
comparison) is shown in Fig. 12. The maximum values in the 
HSS are close to, although not necessarily associated with, 
the climatological probability values indicated for the 
three different time periods. 

Another point needs to be mentioned concerning the HSS. 
The climatological prcbability values for recurvatura (Table 
II) at the specified times should be compared to the means 
and variances of the predictands shown in Table XII. All 
three climatological values of recurvatura fall within the 
first standard deviation around the mean in ail three tests. 
With this in mind, one would predict that methods with the 
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o 




PROBABILITY OF RECURVBTURE 



Figurs 12. 



Held ke Skill Scores with the Value of the 
Predictaad from the PSE Equations that is 
Defined tc be Recurvature Varied n.o 
0.95. 



om 0.05 to 



smallest variances would have a higher HS3 due to the 
greater probability of being near the climatological value 
(low signal regime) than in the tails (high signal regime) 
of the preaictand distribution (to be discussed later). 

In Table XIII, the combined P&E method produced the 
highest scores relative to climatology, except at the 36 h 
period with the dependent sen, where the PERS method was 
best. An indication of the clustering of predictands about 
the mean appears in the 54 h independent set. The 
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predictands based only cn EOFS have the highest scores 
relative to climatology at this time, while at both 36h and 
72 h the EOFS were inferior to the other two methods. P&E 
generally showed the highest skill for the independent set, 
with the exception being at 54 h as discussed above. 

The distribution and cumulative frequency of the 
filtered predictands for the three tests are given in Table 
XIV, XV, and XVI. An unbiased distribution would have 10* 

of the observations in each category. This ideal limit is 
nearly attained at the 36 h time period for the dependent 
data set. However, biases become more evident for later 
time periods. The P&E test generally is less biased than 
the PERS test. 3oth P&E and PSRS are less biased than the 
EOFS test, which again is an indication of the clustering of 
the EOF predictands near the mean value. 

TA3LE XIV 

Predictand Distr ibui 

Dis' 



Predictand Distribution and Cumulative Percentage 
irrbution for 36 h. Cumulative Percentaae is in Parens. 
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TAB LE XV 



Predictand Distribution and Cumulative Percentage 
Distribution for 54 h. Cumulative Percentage is in Parens. 
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TABLE XVI 

predictand Distribution and Cumulative Percentage 
Distribution for 72 h. Cumulative Percentage is in Parens. 
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The Brier P-Scores for the three individual 
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also the BPS for the sample climatology and for chance are 
given in Table XVII which indicates that each of the 
three methods is superior to climatology and chance at 
all forecast times. This is encouraging and is also an 
indication that the SCF approach of describing the synoptic 
forcing is a valid approach that merits further 
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TABLE XVII 



F-Score Statistics for 36, 54 and 72 h Dependent and 
Independent Data Sets After Filtering of the Preaictand 
(A Score of 0.0 is Perfect) 
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indicate the 
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as time increases. It is 


believed that 


this 


decrease in 


skill with respect to time 



could be lessened to some degree by the inclusion cf 
prognostic EOF fields. This was indicated by Shaffer (1982) 
as one of the largest drawbacks of the EOF tropical motion 
forecasting procedure. 
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The third and simplest method of skill analysis is the 
percent correctly forecast (PCr) method. Table XVIII lists 
the FCFs obtained for each test for both the dependent and 
independent data sets. This table again illustrates the 
superiority of the P8E test after the initial 36 h time 
period. Notice again for the small independent data sen at 
54 h, the EOF method is superior to both of the other 
methods which is believed to be due to the clustering 
effect. 



TABLE XVIII 

Percent Co?r^ctly Forecast Events from the Regression Method 
and Defanrng Recurvature as a Predictand Exceeding the 
Climatology Probabilities of Recurvature 
(A Score of 100 is Perfect) 
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This simple method of testing skill was included as an 
attempt to compare the three tested procedures against the 
models currently being used by J TWC Guam for guidance. The 
discussion and comments pertaining to the results of the 



model comparisons are included in Appendix A. 
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VI. CONCLUSION S AND JUSTIFICATION FOR CONTIN UED 

RESEARCH 

Elsberry (1979) discusses the inherent problems and 
advantages of dynamic tropical cyclone models. He states 
that the lack of vertical and horizontal resolution in the 
models, as well as a poor understanding of the environmental 
processes that occur in the vicinity of a tropical storm, 
are limiting the progress toward more accurate forecasts. 
Although nested grid models (Harrison, 1981; Harrison and 
Fiorino, 1982) have been used to improve the resolution, 
there continue to be problems of data acquisition on the 
scales of interest. This problem has no solution in the 

immediate or extended time frames, which necessitates the 
need for other reliable and less costly methods. 

The application of a linear regression scheme to a 
statistical model using a combination of persistence-related 
parameters and an EOF representation of the synoptic forcing 
was completed. The combination of these two types of 
parameters have shown to be superior for forecasting 
tropical cyclone recurvature as defined. 
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The main attributes of the EOF approach to forecasting 

tropical cyclone recurvature and motion (Shaffer, 1982) must 
# 

again be stressed: 

• The method is extremely inexpensive; 

• The forecasts can be done on a handheld programmable 
calculator ; 

• The method has been shown to be more accurate for 
forecasting recurvature than persistence and synoptic 
parameters alone; and 

• The method has shown to give superior forecasts when 
compared to the official JTWC 'forecasts using the 
analysis of the 500 mb fields (Shaffer, 1982). 

It is believed that the inclusion of prognostic synoptic 

fields, as suggested by Shaffer (1982), would improve the 

percentage of the explained variance of the predictand. 

This would also be reflected as giving better results for 

the recurvature study described here. 

An independent operational evaluation of the schemes 

presented here needs to be dons to determine the validity of 

the method. 

The definition of recurvature needs to be modified. A 

tropical storm moving 314° is considered as a 

non- recurvature case where a storm moving 316° is considered 
a recurvature case. Probabilities of re curvature based on 
the departure of the tropical cyclone from a forecast track 
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or a straight line could be utilized instead of using a 
binary definition as done in this research. 

It is understood that the application of new and 
superior methods of forecasting tropical cyclone recurvature 
and motion should be guestionei. It is believed chat this 
change of methodology (ie, using an EOF representation of 
the synoptic fields), if adopted, will improve the image of 
the military weather services by being a positive pioneering 
effort in an area that has been stagnant for many years. 
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APPENDIX A 



MODEL COMPARISONS 



A. EORPOSE 

The purpose of this section is to supply the Western 
Pacific tropical cyclone forecaster with information on the 
various guidance techniques for forecasting tropical cyclone 
motion, and, in particular, their reliability in forecasting 
recurvature. The definition of recur vat ure (Leftwich, 
1978,1982) is similar as described earlier except for the 
northeastward motion during the 24 h (vice 12) prior to 
the end of the forecast period. This modification was done 
to accomodate the' forecast warning verification times which 
are 24, 48 and 72 h. 

"Recurvature is a net disDiacement northward of 315° 
during the forecast oeriod or the attaining o 
northeastward motion during the 24 hours prior to tne en 
of the forecast Deriod. 



B. THE DATA 



The data analyzed 
Naval Environmental 
during the 1979-1981 
data were the best 



were assembled by Dr. Ted Tsui of the 
Prediction Research Facility (NEPRF) 
Northwest Pacific typhoon seasons. The 
track information from the tropical 
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pi Hi 



disturbance stage through the entire life cycle of each of 
the 28 storms of 1979 and 1980 and the 29 storms of 1981. 
Also included were the 24, 48 and 72 h positions of 12 
tropical cyclone forecast aids available at JTWC. Seven of 
these aids and the Joint Typhoon Warning Center official 
forecasts are evaluated below with respect to the 

recurvature decision. 



C. MODEL DESCRIPTIONS COMPARISON TESTS 

The following paragraphs are descriptions of the 
objective forecast aides currently being used at JTWC. The 
^descriptions have been taken from the 1931 Annual Tropical 
Cyclone Report. All JTWC objective techniques, except those 
involving dynamic models, are executed using operational 
data available prior to warning times ( 00GMT and 12GMT). 
These dynamically-based objective aids are usually received 
at JTWC within 2 to 5 hours of a specific warning time. 

An analogue method, abbreviated T0P5, provides three 
movement and intensity forecasts. These three forecasts are 
for analogues of straight-moving cyclones, of recurving 
cyclones and a combination of analogues including straight, 
recurving, and all other cyclones that do not meet the 
specific cirteria of those two categories. 
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An extrapolation technique, abbreviated XTRA, provides a 
track connecting the 12-h old preliminary best track 
position and the current warning position, which is then 
extrapolated to 24 and 43 h. HPAC, which is a combination 
of XTRA and a climatological aid, provides a forecast of the 
track and intensity changes based on the initial position of 
the system. The climatological data base includes the 
period from 1945 to 1973. The ( BPAC) forecasts are based on 
blending the past motion of the tropical cyclone with the 
climatology forecast positions. This program requires only 
a Texas Instrument (TI-59) calculator to generate 12 to 72 h 
forecast positions. The blending routine gives less weight 
to persistence at each succeeding forecast interval. 

An updated version of the HATR ACK/MOHATT steering 

program (now called CYCLOPS and designated as CY50 for the 
500 mb level) can provide steering forecasts at many 

different levels. The program can be run in the modified 
(includes a 12-h persistence bias) or unmodified versions 
applied to either analysis or prognostic fields. The 
program advects a point vortex on a pre-selected analysis 
and/or smoothed prognostic field at designated levels in 6-h 
time steps through 72 h. In the modified version, the 
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program uses the previous 12-h history position to compute 
the 12-h forecast error and applies a bias correction to the 
forecast positions. In 1981, only the modified version, in 
the prognostic mode for the 500 mb level was verified. 

.The dynamic Tropical Cyclone Model (ICMO) is a coarse 
mesh (220 km) primitive equation model. The digitized 
cyclone warning position is bogused into the 350 mb wind and 
temperature fields of the F NOC global band analysis. 
Hemispheric mcdel forecast data are used on the boundaries. 
Two versions are currently run: One runs from forecast 

fields and is available on an automatic request basis, and 
the other is initialized on the analysis fields noted 
earlier. The model based on the analysis fields was tested 
in this study. 

The NTCM or "nested" primitive equation tropical cyclone 
model is initialized on FNOC 1 2-h forecast fields. The 
model covers a limited, but relocatable, tropical domain 
with three layers in the vertical. The finer scale, or 
"nested" grid, covers a 1200 km square area with a 41 km 
grid spacing and moves to keep a 850 mb vortex in its 
center. This model is integrated with a coarse channel 
model grid with a grid spacing of 205 km over a 6400 by 4700 
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km domain. Once initialized, the model runs independent of 
the remaining FNOC forecast fields. The NTCM is available 
via automatic request for 00 and 12 GMT forecast fields 
only . 

D. SC BE ENING PROCEDU BE FOR MODEL TESTING 

Several different requirements were taken into 

consideration for inclusion in this comparison study that 
were not considered for the EOF analysis scheme. Inclusion 
was based on: 

• Using the 'best track' past zonal and meridional 

displacements (after the initial 12 h of the track) 
vice the past warning displacements used in the aOF 
study, to test for eastward displacement (net 

considered case) ; 

• Pseudo- ir. depende nee is again assumed by requiring time 
separations* o f 36 h within the same storm; 

• A 30 h past history was not required; 

• All tropical cyclones were included, rather than iust 
tropical storms and typhoons; and 

• All storms north of the equator, vice 10° north, were 
included . 

Table XIX lists the 24, 48 and 72 h recurvature 

probabilities for the three year sample. when compared to 
the corresponding values of the sample climatology listed in 
Table II, the difference between the inclusion criteria 
become evident. The inclusion of the tropical depressions 
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in the area south of 10° north is the largest contributor to 
the lower values of recurvature probability. A high annual 
variability is indicated in this small sample. There is 
generally a slow increase of probability of recurvature with 
time as was the case in Table II . 

TABLE XIX 



Numbep <?f Recurving (rev) , N on-recurvina (nrev), and 
Probabilxty of Recurvature (prob) for Trbpxcal Cyclone 
Seasons 1979, 1980 and 1981. 







1979 






19 80 






1981 




ta u 


rev 


nrev 


prob 


rev 


nr cv 


prob 


rev 


nrev 


prob 


24 


39 


94 


.293 


36 


7 8 


.316 


26 


73 


.263 


48 


41 


84 


.328 


35 


5 5 


.350 


23 


63 


.267 


72 


25 


81 


.236 


29 


5 3 


.354 


21 


52 


.288 



The HSS for the three separate years are listed in Table 
XX. The HSS in this data set remains as described earlier 
with a slight modification. This difference is tnat the 
verification of the model's predicted recurvature is based 
entirely upon the predictands matching the actual 

recurvature criteria as defined. There is r.o climatological 
threshhold that is set for inclusion as a recurvature or 
non-recurvature as was the case for the earlier data sample. 

Again, the most noticable feature of the table is the 

general decrease of skill with respect to time. The 
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exception to this generality is the 3PAC model performance 
for 1979 and the TCMO model for 1981 where the H5S increased 
with time. An interesting comparison can be shown when 
comparing JTWC with the TCMO forecast aide for 1981. There 
is a definite correlation between TCMO and JTWC. This is an 
indication of JTWC following the guidance of this dynamic 
model quire closely in the extended time frames. 

Table XXI lists the BPS for the three time periods of 
interest. The BPS forecast predictands in this case are not 
on the interval between 1.0 and 0.0, but are identically 
equal to 1.0 or 0.0. Any knowledge of a comparison of BPS 
results using a continuous, predictand as compared to a 
binary predictand scheme is available. Comparison of -he 
common time period of 72 h between the -wo entirely 
independent data sets of different sizes there fore has 
little meaning. The BPS of the TCMO model in 1981 and BP AC 
for 1979 are the only models that approach the BPS for the 
combination of EOFS and Persistence in the earlier study. A 
general trend toward larger BPS are prevalent for all years 
for most of the forecast aids. 

Comparison of the PCF values indicated in Table XXII, 
and the HSS reveals the similarity of the two skill scores. 
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TABLE XX 



Heidke Skill Sco^e for Models for 1979, 1980, and 1931 

(A Score Of 1.0 is Perfect, and a (*) Indicates too 

Few Cases to be Hanked) 



Model 


results for 


19 79 














HP AC 


TCMO 


JTWC 


T0P5 


CY50 


BPAC 


XTRA 


NICM 


24 h 


.449 


.621 


. 424 


.2 96 


.388 


.403 


.304 


.167 


cases 


421 


95 


427 


247 


304 


345 


425 


12 


rank 


(2) 


(D 


(3) 


(7) 


(5) 


(4) 


(6) 


<*) 


48 h 


.457 


. 506 


. 355 


.2 30 


.426 


. 568 


.235 


.667 


cases 


368 


8 1 


363 


2 28 


265 


296 


371 


12 


rank 


(3) 


(2) 


(5) 


(7) 


(4) 


(1) 


(6) 


(*) 


72 h 


.381 


. 391 


. 347 


.266 


.095 


. 586 


.072 


.667 


cases 


252 


69 


288 


1 99 


210 


251 


3 04 


12 


rank 


(3) 


(2) 


(4) 


(5) 


(6) 


(1) 


(7) 


<*) 


Model 


r esul 


ts for 


19 80 














HP AC 


TCMO 


JTWC 


TOPS 


CY50 


BPAC 


XTRA 


NTCM 


24 h 


.5 18 


.621 


. 424 


.296 


.3 88 


.403 


. 304 


.167 


cases 


4 19 


73 


412 


353 


364 


36 1 


373 


62 


rank 


(2) 


(1) 


(3) 


(7) 


(5) 


(4) 


(6) 


(8) 


48 h 


.3 13 


.406 


. 418 


.461 


.4 11 


. 453 


.301 


.387 


cases 


352 


64 


347 


3 04 


309 


300 


306 


62 


rank 


(7) 


(5) 


(3) 


(D 


(4) 


(2) 


(8) 


(6) 


72 h 


.004 


. 378 


. C89 


. 1 12 


.273 


. 2 17 


. 164 


.309 


cases 


261 


45 


257 


245 


245 


240 


225 


55 


rank 


(8) 


(1) 


(7) 


(6) 


(3) 


(4) 


(5) 


(2) 


Model 


r esul 


ts for 


19 81 














HP AC 


TCMO 


JTWC 


TO P5 


CY 50 


BPAC 


XTRA 


CiTCM 


24 h 


.462 


.457 


. 542 


.466 


.567 


. 480 


.433 


.333 


cases 


383 


129 


384 


348 


282 


377 


383 


156 


rank 


(6) 


(7) 


(2) 


(5) 


(D 


(4) 


(3) 


(3) 


48 h 


.298 


. 574 


. 515 


.333 


.439 


. 388 


. 255 


.289 


cases 


325 


108 


309 


297 


235 


320 


3 25 


1 35 


rank 


(6) 


(1) 


(2) 


(5) 


(3) 


(4) 


(8) 


(7) 


72 h 


.0 13 


.647 


. 540 


.021 


.304 


.063 


. 149 


.196 


cases 


235 


68 


268 


235 


1 8 1 


254 


235 


107 


rank 


(8) 


(D 


(2) 


(7) 


(3) 


(6) 


(5) 


(4) 


The rankings 


of these two 


method s 


are 


the same 


for 5 1 


of 6 ' 



comparisons. 
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TABLE XXI 



Bri i57§,' 

a 


Score S 
1980 and 
(*) Indi 


taxistics for 24. 48 
1981 (A Scoca of 0.0 
cates too Few Cases t 


and 72 & Models 
is Perrect, and 
o be Ranked) 


for 


Model 


results for 


19 79 














HP AC 


TCMO 


JTWC 


T0P5 


CY50 


B PA C 


XTRA 


NTCM 


24 h 


.570 


.379 


. 576 


.704 


.612 


.597 


.6 96 


.333 


cases 


4 21 


95 


427 


247 


3 04 


345 


4 25 


12 


rank 


(2) 


(D 


(3) 


(7) 


(5) 


(4) 


(6) 


(*) 


48 h 


.543 


.494 


. 645 


.7 02 


.574 


.432 


.765 


.333 


cases 


368 


81 


363 


2 23 


265 


296 


371 


12 


rank 


(3) 


(2) 


(5) 


(6) 


(4) 


(D 


(7) 


(*) 


72 h 


.6 19 


.609 


. 653 


.734 


.905 


.414 


. 928 


.333 


cases 


252 


69 


288 


1 99 


210 


251 


304 


12 


rank 


(3) 


(2) 


(4) 


(5) 


(6) 


(1) 


(7) 


(*) 


Model 


results for 


19 80 














HP AC 


TCMO 


JTWC 


T0P5 


CY 50 


B P AC 


XTRA 


NTCM 


24 h 


.482 


.493 


. 461 


.465 


.440 


.454 


. 536 


.581 


cases 


4 19 


73 


412 


353 


364 


36 1 


373 


62 


rank 


(5) 


(6) 


(3) 


(4) 


(1) 


(2) 


(7) 


(8) 


48 h 


.688 


. 594 


. 582 


.539 


.589 


.547 


.699 


.613 


cases 


352 


64 


347 


3 04 


309 


300 


306 


62 


rank 


(7) 


(5) 


(3) 


(D 


(4) 


(2) 


(8) 


(6) 


72 h 1 


.004 


.622 


. 911 


.873 


.727 


. 783 


. 836 


.691 


cas es 


261 


45 


257 


2 45 


245 


240 


225 


55 


rank 


(8) 


(1) 


(7) 


(6) 


(3) 


(4) 


(5) 


(2) 


Model 


results for 


1981 














HP AC 


TCMO 


JTWC 


T0P5 


CY 50 


3PAC 


XTRA 


NTCM 


24 h 


.538 


.543 


. 458 


.534 


.433 


. 520 


.517 


.66 7 


cases 


383 


129 


384 


348 


282 


377 


383 


156 


rank 


(6) 


(7) 


(2) 


(5) 


(1) 


(4) 


(3) 


(8) 


48 h 


.702 


.426 


. 485 


.667 


.511 


.612 


.745 


.711 


cases 


325 


108 


309 


297 


235 


320 


3 25 


135 


rank 


(6) 


(1) 


(2) 


(5) 


(3) 


(4) 


(8) 


(7) 


72 h 


.9 87 


.35 3 


. 460 


.9 79 


.6 96 


.937 


.902 


.804 


cases 


235 


68 


26 8 


235 


181 


254 


235 


1 07 


rank 


(8) 


(D 


(2) 


(7) 


(3) 


(8) 


(5) 


(4) 
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TABLE XXII 



Percent Correct Forecast Statistics for Models for 1979, 
1980 and 1981 (A Score Of 100 is Perfect t and a (*) 
Indicates t co Few Cases to be Ranked) 



Model results for 1979 





HP AC 


TCMO 


JTWC 


T0P5 


CY 50 


B P AC 


XTRA 


NTCM 


24 h 


71 .5 


81. 1 


7 1.2 


64.8 


69.4 


70. 1 


65.2 


58.3 


cases 


421 


95 


427 


2 47 


3 04 


345 


425 


12 


rank 


(2) 


0) 


(3) 


(7) 


(5) 


(4) 


(6) 


(*) 


48 h 


72.8 


75. 3 


6 7.8 


64.9 


71 .3 


78.4 


61.7 


83.3 


cases 


368 


81 


363 


2 28 


255 


296 


371 


12 


rank 


(3) 


(2) 


(5) 


(6) 


m 


(1) 


(7) 


(*) 


72 h 


69.0 


69. 6 


6 7.4 


63.3 


54 .8 


79.3 


53.6 


83.3 


cases 


252 


69 


283 


1 99 


210 


251 


304 


12 


rank 


(3) 


(2) 


<*) 


<5) 


(6) 


(D 


(7) 


(*) 


Model 


. results for 
HP AC TCMO 


19 80 
JTWC 


TOP5 


CY 50 


B P AC 


XTRA 


NTCM 

71. 0 


24 h 


75 . 9 


75.3 


7 6.9 


76.8 


78.0 


77.3 


73.2 


cases 


419 


73 


412 


353 


364 


361 


373 


62 


ran k 


(5) 


(6) 


(3) 


(4) 


(D 


(2) 


(7) 


(8) 


48 h 


65.6 


70. 3 


7 0. 9 


64.9 


70 .6 


72. 7 


69.3 


69.4 


cases 


352 


64 


347 


304 


309 


300 


306 


62 


rank 


(7) 


(4) 


(2) 


(8) 


(3) 


(1) 


(6) 


(5) 


72 h 


49.8 


68.9 


5 4. 5 


63.3 


63 .7 


6 0.8 


64.0 


65.4 


case s 


26 1 


45 


257 


245 


245 


240 


225 


55 


rank 


(3) 


(1) 


(7) 


(5) 


(4) 


(6) 


(3) 


(2) 


Model 


. result 
HP AC 


.s for 
TCMO 


19 81 
JTWC 


T0P5 


CY 50 


3 ? AC 


XTRA 


NTCM 


24 h 


73. 1 


72. 9 


7 7. 1 


73.3 


78.4 


74.0 


74.2 


66.7 


cases 


383 


129 


384 


348 


282 


377 


383 


156 


rank 


(6) 


(7) 


(2) 


(5) 


(1) 


(4) 


(3) 


(8) 


48 h 


64 .9 


78.7 


7 5.7 


66.7 


74 . 5 


69.4 


62.3 


64 . 4 


cases 


325 


103 


309 


297 


235 


320 


3 25 


135 


rank 


(6) 


(D 


(2) 


(5) 


(3) 


(4) 


(8) 


( 7 ) 


72 h 


50.6 


83.4 


7 7.0 


51 . 1 


65. 2 


53. 1 


54.9 


59.3 


cases 


235 


68 


268 


235 


181 


254 


235 


107 


rank 


(3) 


(D 


(2) 


(7) 


(3) 


(6) 


(5) 


(4) 


E. 


SUMMA3 

Seven 


Y 

for 


ecast 


aids 


used 


by J o i 


nr 


Ty phoc 



Warning Center (JTWC) Guam, were tested for ability to 
forecast tropical cyclone recurvature at 24, 48 and 72 h. 
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Using the skill scores described in Chapter V for the 
1979-81 typhoon season, model performance indicates a large 
annual variation of skill between the models. 
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